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Bispectriun  of  ship-rudiated  noise 


B.J.  Sullivan,  M.J.  Hinich  and  D,  Marandino 


Executive  Summary!  Passive  sonar  performance  is  limited  by  back¬ 
ground  noise  masking  the  noise  radiated  by  the  target  vessel  because  both 
kinds  of  noise  can  appear  similar  to  present  signal- processing  techniques. 
Even  in  the  case  where  the  target  can  be  identified  by  characteristic  narrow- 
band  specttal  lines,  a  sufficiently  high  signal- to* noise  level  is  necessary  for 
detection.  This  report  presents  a  different  approach,  known  as  bispectral 
analysis,  which  provides  a  potential  improvemnet  in  the  ability  to  discrim¬ 
inate  target  noise  from  background  noise  by  utilizing  a  hidden  feature  of 
ship-generated  noise  that  is  not  present  in  the  background  noise. 

Due  to  nonlinear  mechanisms  of  ship  noise  generation,  there  are  inteimod- 
ulations  between  the  different  frequency  components  of  the  noise.  These 
intermodulations  can  be  identified  by  bispectral  analysis.  Instead  of  having 
a  single  frequency  coordinate  as  in  the  conventional  spectrum,  the  bispec¬ 
trum  has  two.  The  bispectrum  therefore  is  represented  by  a  surface  con¬ 
tour,  whereas  the  conventional  spectrum  is  represented  by  a  one-dimensional 
form.  It  is  the  shape  of  this  surface  that  indicates  the  presence  of  target¬ 
generated  noise  through  its  intermodulations.  Any  significant  deviation  from 
flatness  of  this  surface  signifies  the  presence  of  nonlinear  components  in  the 
noise-generating  mechanism.  The  hispectrum  therefore  has  the  potential  of 
providing  a  large  signal-processing  gain  in  recognizing  shipping  noise  in  a 
strong  ambient  background,  which  could  provide  a  significant  increase  in 
the  sonar  detection  range. 

The  results  of  a  bispectrum  calculation  of  the  radiated  noise  of  a  ship  re¬ 
ceived  on  a  towed  array  being  towed  by  the  same  ship  is  presented  in  this 
report.  The  array  was  beamformed  so  that  the  forward  endflze  beam  pointed 
at  the  ship  and  the  broadside  besm  sampled  the  environment.  In  this  man¬ 
ner  the  ship  could  effectively  be  ‘switched’  on  and  off.  The  results  show 
that  the  ship  noise  had  significantly  higher  bispectral  components  than  the 
ambient  noise. 

Although  this  work  verifies  the  identification  of  the  ship  through  its  non¬ 
linear  features,  the  details  of  the  relationship  between  the  bispectrum  itself 
and  the  particular  nonlinear  mechanisms  remains  to  be  investigated. 
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Bispectrum  of  ship-radiated  noise 
E.J.  Sullivan,  M.J.  Hinich  and  D,  Marandino 


Abstracts  The  bispectrum  of  ship-radiated  noise  is  estimated.  The  noise 
was  received  on  a  towed  array  being  towed  by  the  same  ship  that  served  as 
the  noise  source.  The  array  was  beamformed  such  that  the  forward  endfire 
beam  pointed  toward  the  towing  platform  and  the  broadside  beam  sampled 
the  environment  without  the  radiated  noise  of  the  ship.  The  results  show 
that  there  exist  frequency-dependent  bispectral  components  in  the  ship’s  ra¬ 
diated  noise,  whereas  the  ambient  noise  does  not  contain  any  significant  bis¬ 
pectral  components.  Since  the  existence  of  a  non-sero  frequency-dependent 
bispectrum  indicates  the  existence  of  nonlinear  components  in  the  noise- 
generating  mechanism,  it  is  concluded  that  the  radiated  noise  of  the  towing 
platform  contains  such  nonlinear  mechanisms.  Therefore,  the  bispectrum 
could  be  used  to  indicate  the  existence  of  such  noise  sources  as  would  nor¬ 
mally  be  hidden  in  the  background  noise  when  the  usucj  spectral  estimation 
procedures  are  applied. 

Keywords:  ambient  noise  o  bispectrum  o  gaussianity  o  Hinich  test 
o  nonlinearities  o  ship- radiated  noise  o  transients 
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1.  Introduction 

Recent  work  which  uses  bispectral  analysis  to  characterize  ocean  acoustic  noise  gives 
strong  evidence  that  shipping-dominated  noise  contains  strong  nonlinear  components 
in  its  generating  mechanisms  whereas  the  ambient  sources  do  not  (1).  The  value  of 
the  bispectrum  is  that  it  is  able  to  distinguish  via  the  Hinich  test  [2]  whether  or 
not  the  source  of  the  data  contains  such  nonlinear  components.  Standard  statistical 
tests  are  unable  to  make  such  a  distinction.  Since  these  results  suggest  a  means 
for  differentiating  between  shipping  noise  and  at  least  sc.ne  of  the  other  forms  of 
ambient  noise  sources,  it  is  of  interest  to  make  a  positive  identification  between  a 
specific  source  of  ship-generated  noise  and  its  bispectrum.  Such  an  identification  is 
made  in  this  work  where  we  analyze  data  taken  from  a  towed  array  where  its  towing 
platform  served  as  the  noise  source.  The  array  was  beamformed  in  such  a  manner 
that  it  was  possible  to  'switch’  the  source  on  and  off.  That  is,  bispectra  were  made 
of  the  records  obtained  with  the  main  beam  steered  to  forward  endfire  (toward  the 
towing  platform)  and  compared  to  bispectra  of  the  records  obtained  when  the  main 
beam  was  at  broadside.  The  bispectrum  of  the  forward  endfire  data  indicated  a 
significantly  larger  content  of  nonlinearity  than  that  of  the  broadside  data. 
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2.  Nonlinear  mechanisms  and  models 


A  time  series  is  linear  if  it  can  be  modeled  in  the  form 

x(t)  =  /i(m)<(f  -  m),  (1) 

m 

where  e(n)  is  a  purely  random  sequence,  i.e,  it  is  made  up  of  mutually  independent 
elements.  Thus,  a  white  gaussian  process  is  purely  random  as  a  consequence  of  the 
fact  that  it  is  uncorrelated  and  independent,  If  the  process  is  not  gaussian,  the 
whiteness  is  not  a  sufficient  condition  for  it  to  be  purely  random.  A  generalization 
of  Eq.  (1),  which  forms  a  class  of  nonlinear  models,  is  the  Volterra  expansion  [3], 
which  is  given  by 

»(0 = hm ~ ,n)  +  XI  *•».»«(*  -  m)((i  ~ n) 

m  m  n 

+  £  S  Am,»lPe(f  -  -  n)f(<  -  p) 

m  n  p 

+  ■  ■  •  (2) 

where  the  hm  n.  are  called  the  Volterra  kernels,  The  higher-order  terms  in  Eq.  (1) 
represent  the  nonlinearities  which  would  not  be  explicitly  manifested  by  standard 
tests.  An  explicit  example  of  this  is  the  bilinear  model  considered  by  Priestly  [4] 
which  is  given  by 

x(t)  =  e(t)  +  ae(t-  2)<(t  -  1),  (3) 

where  c  is  purely  random.  It  is  easy  to  show  that  the  variance  of  x(t)  is  given  by  the 
variance  of  t  and  x(t)  is  uncorrelated.  Therefore  any  whiteness  test  on  x(t)  would 
confirm  that  it  is  white  and  the  erroneous  conclusion  would  be  drawn  that  there  ic 
no  remaining  model  structure  to  fit. 

From  a  heuristic  point  of  view,  since  the  major  effect  of  the  nonlinearities  is  to 
cause  intermodulation  between  the  frequency  components  of  the  driving  process, 
without  phase  information  in  addition  to  the  spectral  power,  the  presence  of  the 
nonlinearities  will  not  be  detected. 

It  can  be  seen  from  Eq.  (2)  that,  even  for  relatively  simple  cases,  there  can  be  a 
large  number  of  parameters  to  deal  with  in  studies  of  nonlinear  models.  For  this 
reason,  most  studies  treat  special  cases  as  nonlinear  AR,  MA  or  ARMA  models  as 
in  the  example  given  above.  The  bispectra  of  a  representative  set  of  these  models 
were  computed  by  Ashley  et  al.  (5|. 
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Generally  speaking,  ship-radiated  noise  would  be  expected  to  arise  from  two  sources: 
hydrodynamic  and  mechanical.  The  major  hydrodynamic  source  of  ship-radiated 
noise  is  propeller  cavitation,  which  can  be  quite  broadband  and  falls  off  roughly  as 
the  square  of  the  frequency  [6j.  Cavitation  noise  by  its  very  nature  is  a  nonlinear 
phenomenon  since  it  arises  from  the  collapse  of  bubbles  which  can  happen  in  a 
very  nonlinear  way.  Direct  radiation  from  other  hydrodynamic  sources,  such  as  flow 
noise,  would  not  be  expected  to  make  a  strong  contribution  to  the  energy  in  the  sonar 
frequency  range.  However,  flow  will  play  a  major  role  by  coupling  energy  into  the 
ship’s  hull  which  then  can  radiate  acoustic  energy  into  the  water.  Even  in  the  case  of 
small  amplitude  vibration  of  hull  plates,  any  structural  characteristics  that  produce 
a  nonlinear  Hooke’s  law  phenomenon  would  constitute  a  nonlinear  mechanism.  From 
this  point  of  view  then  it  would  be  expected  that  along  with  cavitation  the  nonlinear 
characteristics  of  the  mechanical  sources  of  noise  such  as  hull  plates  and  rotating 
machinery  would  also  play  a  role  in  the  bispectrum  of  ship-radiated  noise. 

There  do  not  seem  to  be  a  great  deal  of  explicit  examples  in  the  recent  literature 
connecting  such  physical  processes  to  nonlinear  time  series.  In  an  attempt  to  provide 
some  motivation  for  the  analysis  contained  in  this  work  we  point  out  that  some  of  the 
models  in  [5]  can  be  identified  with  mechanical  systems.  In  particular,  the  example 
referred  to  as  the  ‘threshold  AR’  in  [6]  given  by 

*(<)  -/*»«  -  1) +  €(<). 

where 

-0.5  if  *(<  -  1)  <  1 
0.4  otherwise, 

can  be  identified  with  an  approximation  to  an  overdamped  harmonic  oscillator  with 
a  threshold  value  for  its  damping  coefficient.  Another  example  that  identifies  non¬ 
zero  biapectra  with  mechanical  sources  is  given  by  Sato  et  al.  [7],  where  gear  noise  is 
analyzed.  Although  they  do  not  give  a  specific  example  of  a  time  series  associated 
with  their  physical  model,  the  efficacy  of  bispectral  analysis  in  identifying  such 
mechanical  sources  of  acoustic  noise  is  demonstrated. 
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3.  The  bispectrum  as  a  test  for  nonlinearity  and  nongaussianity 

The  bispectrum  of  a  time  series  i(t)  is  defined  to  be  the  double  fourier  transform  of 
the  third-order  cumulant,  i.e. 

»(«,.«,)  =  EH(7*-(m.n)e-<<u'im+"in),  (5) 

m  n 

where 

C-mi  =  T  m)x(t  +  n)*(t)).  (6) 

If  r.(t)  can  be  represented  by  Eq.  (1)  it  follows  that  its  power  spectrum  is  given  by 

S(ui)  =  crJ|G(ui)|1  (7) 

with 

G(u>)  =FT{*(<)}  (8) 

<T*  =£?(«*].  (9) 

Thus  if  sc(<)  is  linear,  the  bispectrum  can  be  expressed  in  the  form 

B(ui ,  uij)  =  pjGJuq  )G(u>j  )G'(a>i  +u>j),  (10) 

where  is  the  third  moment  of  t,  i.e. 

=  £(*»].  (11) 

The  bispectrum  has  a  principal  domain  that  is  triangular  This  follows  from  its 
symmetries.  Letting  fj  and  fk  label  the  two  frequencies,  the  principal  domain  is 
defined  by  the  inequalities  0  <  j  <  j-lV  and  0  <  k  <  j,  where  2 j  +  k  <  N  with  N 

being  the  size  of  the  data  set.  This  principal  domain  is  shown  in  Fig.  1  where  it 

can  be  seen  that  it  is  divided  into  two  regions  which  are  labelled  IT  (inner  triangle) 
and  OT  (outer  triangle).  Hinich  has  shown  that  if  the  sample  rate  of  the  data  is 
equal  to  or  exceeds  the  Nyquist  rate,  then  the  bispectrum  in  the  OT  is  expected  to 
be  zero  [8].  The  appearance  of  non  zero  values  of  the  bispectrum  in  the  OT  is  an 
indication  of  aliasing,  and  therefore  transient  s. 

The  Hinich  test  is  based  on  a  prewhitened  bispectral  gain  called  the  nkcwntss  func¬ 
tion  which  is  given  bv 

(  '’  j)  (5(ui,)5(uiI)5(u»I  +«,)]*/»•  (  ) 
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Fig.  i.  The  frequency  plane  lattice  in  the  principal 
domain.  N  is  the  block  length. 


From  Eq.  (7)  and  (10),  we  see  that  if  x(t)  is  linear  it  immediately  follows  that 

r*(uq,W2)  =  nl/crS.  (13) 

For  a  linear  process  then,  the  skewness  is  constant  over  t.,e  frequency  plane.  Note 
that  since  fi3  —  0  for  the  gsussian  case,  a  non-zero  value  of  the  bispectrum  rejects 
gaussianity.  Consequently,  all  nonlinear  processes  axe  non-gaussiap. 

In  order  to  test  for  nongaussianity,  Hinich  compares  a  normalized  estimate  of 
T  (uq,w2)  at  different  points  on  the  w-plane  to  a  central  chi  squared  distribution 
with  two  degrees  of  freedom.  The  chi-squared  statistics  are  summed  over  a  grid  of 
bifrequencies  (frequency  pairs)  in  the  the  principal  domain.  Any  significant  devia¬ 
tion  of  the  sum  of  these  statistics  from  a  chi-square  of  2 P  degrees  of  freedom  is  taken 
to  be  a  rejection  of  the  gaussian  assumption,  where  P  is  the  number  of  bifrequencies 
in  the  principal  domain, 

The  test  used  by  Hinich  for  linearity  is  based  on  the  dispersion  of  the  estimate 
of  the  chi-squared  statistics  assuming  a  constant  noncentrality  parameter  which  is 
estimated  from  the  data.  The  noncentrality  parameter  for  each  chi  squared  statistic 
is  the  squared  skewness  function  at  that  bifrequency  scaled  by  the  normalization 
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constant.  Thus  if  the  process  is  Linear,  the  skewness  is  constant,  implying  that  all 
the  noncentrality  parameters  are  equal.  This  dispersion  is  computed  for  a  given 
fractile  range  and  compared  to  the  dispersion  over  the  same  fractile  range  for  a 
gaussian  process.  If  this  dispersion  significantly  exceeds  that  of  the  linear  process, 
the  the  assumption  of  linearity  is  rejected.  Hinich  originally  used  the  inner  quartile 
range  for  this  test*.  However,  here  we  shall  use  the  80th  percentile  range  since  it 
seems  to  be  quite  robust.  For  details  of  these  tests  see  the  appendix  in  [1). 
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4.  Data  characteristics 

As  mentioned  in  the  introduction,  the  data  were  taken  from  a  towed  array.  The  tow¬ 
ing  platform  was  the  Mana  Paoltna  G.,  the  former  research  vessel  of  the  SACLANT 
Undersea  Research  Centre.  It  has  a  length  of  79  m,  a  beam  of  12  tn  and  a  draft 
of  4.5  m.  Its  displacement  is  2800  metric  tonnes.  The  engine  is  a  six  cylinder  two- 
stroke  single  action  diesel  with  1800  BHP.  At  a  speed  of  5  kn  the  shaft  speed  is 
110-115  turns  per  minute.  The  fixed-pitch  propeller  has  5  blades. 

The  array  was  towed  at  a  speed  of  5  kn  at  a  depth  of  70  m,  with  a  scope  of  750  ill. 
This  placed  it  at  the  top  of  the  thermocline.  There  was  a  weak  sound  channel  at 
100  m.  The  area  was  the  Ionian  basin  region  of  the  Mediterranean  where  the  depth 
was  about  2500  m.  The  weather  was  good  with  a  sea-state  of  2.  There  was  no 
shipping  in  the  immediate  area 

The  data  consist  of  two  records,  each  taken  at  different  times  during  the  same  day, 
which  will  be  referred  to  as  File  1  and  File  2.  Each  record  contains  50  s  of  endfire  data 
and  50  s  of  broadside  data.  The  spectrum  had  a  flat  region  about  120  Hz  wide  with 
a  center  frequency  of  1090  Hz.  These  spectrum  characteristics  are  a  consequence 
of  the  fact  that  the  data  were  filtered  with  a  linear  phase  Parks-McClellan  filter 
with  127  coefficients  and  an  out-of-band  attenuation  of  54  dB.  The  spectra  for  the 
broadside  and  endfire  cases  for  File  1  are  shown  in  Fig.  2.  The  towed  array  was 
made  up  of  40  elements  with  an  interelement  spacing  of  49  cm.  Thus  the  data  are 
spatially  oversampled. 

The  two  data  records  consist  of  slightly  mor'-  than  30  800  observations.  These  records 
were  divided  into  blocks  and  then  averaged  over  these  blocks.  Two  block  lengths 
were  used  in  the  actual  calculation.  These  were  of  length  140  and  160.  Since  the 
sample  interval  was  0.0016  s,  the  respective  bandwidths  were  4.46  Hz  and  3.85  Hz. 
The  number  of  observations  allowed  an  average  over  220  blocks  for  the  block  length 
of  140  and  an  average  over  192  blocks  for  the  block  of  length  160.  The  spectra  and 
bispectra  were  then  computed  using  these  averaged  blocks. 
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5.  Results 

As  was  pointed  out  in  Sect.  3,  the  bispectrum  has  a  triangular  principal  domain. 
This  is  divided  into  an  inner  triangle  (IT)  and  an  outer  triangle  (OT)  where  the 
bispectrum  is  zero  in  the  outer  triangle  for  the  case  of  Nyquist  sampling.  The 
bispectrum  in  this  domain  is  thus  expected  to  be  non-zero  in  the  case  of  transients. 
In  Table  1,  the  statistic  which  tests  for  gaussianity  is  tabulated,  both  for  the  complete 
principal  domain  and  the  outer  triangle  only.  This  statistic  is  based  on  the  chi- 
squared  distribution  function  for  the  case  of  2 P  degrees-of- freedom,  where  P  is 
the  number  of  frequency  pairs  in  the  domain  under  test.  Thus,  under  the  null 
hypothesis  of  gaussianity,  the  statistic  would  be  expected  to  be  a  sample  from  such 
a  distribution.  The  statistic  is  presented  in  units  of  the  standard  deviation  under  the 
gaussian  assumption.  It  can  be  seen  from  Table  1  that  in  all  cases  the  endfire  data  in 
the  principal  domain  are  strongly  non-gaussian.  For  the  case  of  the  broadside  data 
in  the  principal  domain,  only  the  data  from  File  1  are  strongly  non-gaussian.  In  the 
outer  triangle,  the  endfire  data  from  both  data  files  are  strongly  non-gaussian  and 
the  broadside  data  from  both  files  are  not.  These  conclusions  hold  for  both  block 
sizes. 


ruble  1 

Gaussianity  statistics  for  the  two  block  lengths  of  140  and  161) 


Data  set1 

Principal  domain 

Outer  triangle 

Block  size  =  160 

B1 

3,09 

0.28 

El 

13.15 

7.86 

B2 

1.89 

0.06 

E2 

30.77 

13.50 

Block  site  =  140 

Bl 

4.02 

1.08 

El 

13.29 

5.36 

B2 

0.36 

1.79 

E2 

24.84 

11.84 

1  B1  and  El  correspond,  respectively,  to  the  broadside  and  endfire 
cases  for  File  1.  Similarly.  B2  and  E2  correspond  to  broadside  and 
endfire  for  File  2. 


The  statistic  which  tests  for  nonlinearity  is  the  amount  by  which  the  80th  percentile 


-  9  - 


SaCLANTCEN  SM-201 


exceeds  the  same  percentile  under  the  null  hypothesis  that  the  skewness  is  constant, 
which  is  implied  by  a  linear  hypothesis.  It  is  in  units  of  the  standard  deviation 
under  the  gaussian  assumption.  This  statistic  is  tabulated  in  Table  2.  As  in  the 
case  of  the  gaussian  statistic,  we  present  the  results  for  the  principal  domain  and 
also  for  the  outer  triangle  only.  Here  it  can  be  seen  that  the  endfire  data  from  both 
files  and  for  both  block  sizes  are  strongly  nonlinear.  In  the  case  of  the  broadside 
data,  the  statistic  indicates  a  small  but  statistically  significant  nonlinear  component, 
especially  in  the  case  of  the  data  from  File  1.  This  is  most  likely  due  to  leakage  of 
ship-radiated  noise  into  the  sidelobes  of  the  array  beam  pattern.  As  can  be  seen 
from  Fig,  1,  the  endfire  data  are  slightly  more  than  50  dB  above  the  broadside  data 
in  level.  Therefore,  since  the  sidelobe  rejection  ratio  could  not  be  expected  to  be 
anywhere  near  this  amount,  it  is  not  surprising  to  find  such  leakage.  It  would  be 
difficult  to  estimate  the  actual  rejection  ratio  of  the  array  in  this  situation  since  the 
true  angle  of  the  radiation  from  the  towing  platform  with  respect  to  the  array  was 
not  known. 


Table  2 

Linearity  statistics  for  the  two  block  lengths  of  140  and  160 


Data  set1 

Principal  domain 

Outer  triangle 

Block  size  —  160 

Bl 

2.32 

0.63 

El 

9.68 

5.54 

B2 

1.25 

0.21 

E2 

14.93 

6.21 

Block  size  =  140 

Bl 

2.13 

1.53 

El 

7.63 

3.35 

B2 

-0.22 

-1.84 

E2 

11.14 

6.12 

1  B1  and  El  correspond,  respectively,  to  the  broadside  and  endfire 
cases  for  File  1.  Similarly,  B2  and  E2  correspond  to  broadside  and 
endfire  for  File  2. 


In  order  to  give  some  feeling  for  the  bispectral  levels  associated  with  the  statistics 
presented  in  Tables  1  and  2,  The  bispectra  for  some  selected  cases  along  with  their 
cumulative  distribution  levels  under  the  gaussian  assumption  are  shown  in  Figs.  3 
through  10.  Only  cases  based  on  a  block  size  of  160  are  presented  since  there  was  no 
significant  difference  between  these  results  and  those  for  the  block  size  of  140.  The 
major  differences  occurred  between  the  two  data  files,  which  would  indicate  some 
degree  of  non-stationarity  in  the  data.  Figs.  3  and  4  depict  the  bispectruin  and  the 
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associated  cumulative  distribution  level  for  the  case  of  the  broadside  data  from  File  1 . 
In  Fig.  4  and  all  similar  figures  that  follow,  only  the  cumulative  distribution  level 
above  the  1%  levels  are  shown.  In  Fig,  5  and  6,  the  endfire  data  corresponding  to 
Figs,  3  and  4  are  depicted.  Figures  7  and  8  show  the  bispectrum  and  the  associated 
cumulative  distribution  level  under  the  gaussian  assumption  for  the  case  of  the 
broadside  data  from  File  2,  and  Figs.  9  and  ID  show  the  endfire  case  corresponding 
to  Figs,  7  and  8.  As  can  be  seen  by  comparing  these  figures  with  their  corresponding 
statistics  in  Tables  1  and  2,  it  is  difficult  to  draw  any  substantial  conclusions  directly 
from  the  bispectra  themselves. 
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6.  Conclusions  and  recommendations 

The  major  conclusion  to  be  drawn  from  this  study  is  that  ship-radiated  noise  can  be 
expected  to  contain  a  significant  component  of  noise  emanating  from  sources  that 
contain  nonlinear  characteristics,  and  the  bispectrum  provides  a  means  of  detecting 
such  sources.  The  results  contained  in  this  present  study  can  be  considered  to  be 
particularly  significant  since  the  data  were  reasonably  narrow  band,  with  a  band¬ 
width  of  about  130  Hi  and  a  center  frequency  of  1190  Hi  -  which  corresponds  to 
a  bandwidth-frequency  ratio  of  about  11%.  Thus  a  great  deal  of  the  broadband- 
radiated  noise  was  not  considered  by  the  bispectrai  calculation.  In  particular,  most 
of  the  cavitation  noise,  which  is  expected  to  be  the  major  contribut  or  to  the  radinted 
noise,  will  lie  outside  this  particular  hand. 

A  second  conclusion  is  that  there  appears  to  be  a  significant  amount  of  transients 
in  the  data.  If  it  were  to  be  demonstrated  that  this  is  characteristic  of  ship-radiated 
noise  in  general,  then  this  would  serve  as  another  indicator,  along  with  the  nonlinear 
character  of  the  source,  to  identify  ship-radiated  noise  in  a  strong  background  of 
ambient  acoustic  noise  from  other  sources.  Thus,  the  potential  of  bispectrai  analysis 
in  obtaining  an  effective  increase  in  processing  gain  in  the  detection  of  shipping  noise 
in  a  strong  ambient  background  is  self-evident. 

Another  important  point  brought  out  by  the  results  of  this  study  is  that  in  the 
application  of  the  bispectrum  to  stochastic  data,  such  as  has  been  done  here,  it  is 
important  to  rely  on  well-thought-out  statistical  indicators  rather  than  direct  ob¬ 
servation  of  the  bispectrai  plots  themselves,  as  can  be  seen  by  comparing  Figs,  3 
through  10  to  their  corresponding  statistics  tabulated  in  Tables  1  and  2.  However, 
in  this  regard  it  should  be  added  that  the  configuration  of  the  bispectrai  plots  them¬ 
selves  does  contain  information  as  to  the  particular  type  of  nonlinearities  associated 
with  the  data.  Further  studies  in  this  direction  would  be  of  interest,  particularly  in 
regard  to  the  problem  of  sonar  classification.  In  this  same  vein,  it  would  also  be  of  in¬ 
terest  to  look  at  higher-order  spectra,  since  there  are  certain  classes  of  nonlinearities 
that  do  not  manifest  themselves  in  the  bispectrum  (5|. 
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